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Growth and calcification of the marine coccolithophorid Emiliania huxleyi is affected 
by ocean acidification and macronutrients limitation and its response varies between 
strains. Here we investigated the physiological performance of a highly calcified E. huxleyi 
strain, NZEH, in a multiparametric experiment. Cells were exposed to different CO2 
levels (ranging from 250 to 1314u,atm) under three nutrient conditions [nutrient replete 
(R), nitrate limited (-N), and phosphate limited (-P)]. We focused on calcite and organic 
carbon quotas and on nitrate and phosphate utilization by analyzing the activity of 
nitrate reductase (NRase) and alkaline phosphatase (APase), respectively. Particulate 
inorganic (PIC) and organic (POC) carbon quotas increased with increasing CO2 under 
R conditions but a different pattern was observed under nutrient limitation. The PIC:POC 
ratio decreased with increasing CO2 in nutrient limited cultures. Coccolith length increased 
with CO2 under all nutrient conditions but the coccosphere volume varied depending on 
the nutrient treatment. Maximum APase activity was found at 561 jxatm of CO2 (pH 7.92) 
in -P cultures and in R conditions, NRase activity increased linearly with CO2. These results 
suggest that E. huxleyi's competitive ability for nutrient uptake might be altered in future 
high-C02 oceans. The combined dataset will be useful in model parameterizations of the 
carbon cycle and ocean acidification. 
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INTRODUCTION 

Since the beginning of the industrial revolution, atmospheric 
CO2 has increased at the fastest rate experienced by the Earth 
in the last 65 million years (Zachos et al., 2001). Consequently, 
the increase in CO2 dissolution in seawater has been altering the 
balance of the inorganic carbon species leading to a decrease in 
pH predicted to intensify over the next century (Raven et al, 
2005). Additionally, global warming derived from an increase 
in greenhouse gases induces stratification of the water column 
reducing mixing processes that maintain nutrient supply into the 
euphotic zone (Sarmiento et al., 1998). Therefore, the duration 
and timing of nitrate and phosphate limitation and the bio- 
geographical regions affected are likely to vary in future oceans 
(Sarmiento et al., 1998). Coccolithophores play a major role 
in the carbon cycling being an important source of calcite in 
the open ocean (Gehlen et al, 2007). Emiliania huxleyi is the 
most abundant bloom-forming coccolithophore species (Tyrrell 
and Merico, 2004), and its calcification process is known to be 
affected by variations in carbon chemistry (Riebesell et al, 2000; 
Iglesias-Rodriguez et al., 2008; Langer et al., 2009). However, this 
response varies if other environmental parameters such as nutri- 
ent availability, temperature or light are simultaneously changed 
(Zondervan et al., 2002; Sciandra et al., 2003; Feng et al, 2008; De 
Bodt et al, 2010; Borchard et al, 2011). 



E. huxleyi flourishes after the demise of diatoms, when sili- 
cate, nitrate, and phosphate are limiting (Litchman et al., 2006). 
This ecological strategy stems from a low nutrient quota and 
an extremely high phosphate affinity under phosphate-limiting 
conditions (Riegman et al, 2000). E. huxleyi also takes up nitro- 
gen compounds other than nitrate (Benner and Passow, 2010; 
Bruhn et al, 2010), and assimilates nutrients from organic sources 
through the controlled expression of enzymes active in specific 
metabolic pathways (Dyhrman and Palenik, 2003; Bruhn et al., 

2010) . Despite the sensitivity of E. huxleyi to CO2 (Riebesell et al, 
2000; Iglesias-Rodriguez et al, 2008; Langer et al, 2009) and its 
diversity of nutrient acquisition pathways, the majority of con- 
temporary work has not considered the effect of [CO2] on the 
efficiency of nutrient assimilation. Previous studies have inde- 
pendently assessed the biogeochemical responses of E. huxleyi to 
high CO2 under nitrogen (Sciandra et al., 2003; Leonardos and 
Geider, 2005; Miiller et al., 2012) or phosphorus (Borchard et al., 

2011) limitation using three different E. huxleyi strains (TW1 
PML, B92/11 and a strain isolated in the Raunefjord, Norway). 
It is well accepted that different E. huxleyi strains respond differ- 
ently to varying CO2 levels (e.g., Langer et al, 2009). For example, 
the E. huxleyi NZEH strain presents contrasting calcification 
responses to elevated CO2 compared to other strains tested in the 
laboratory. Similarly, a recent field study revealed the presence 
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of a heavy calcified E. huxleyi morphotype (R-morphotype) in 
"acidic" waters of the South Pacific Ocean, as an exception to 
the global correlation found between coccolithophore calcifica- 
tion and CO^ 2 concentration (Beaufort et al., 2011). Considering 
this biological variability and the projected changes in the extent 
of oligotrophic waters, studying the effect of nutrient limitation in 
conjunction with ocean acidification in different E. huxleyi strains 
is crucial. 

In this study, we assessed the combined effects of elevated 
atmospheric CO2 and nitrate or phosphate limitation on the 
physiology of the E. huxleyi strain NZEH. This is a highly calcified 
strain bearing coccoliths that display the R-morphotype whose 
production appears to be resilient to ocean acidification (Iglesias- 
Rodriguez et al., 2008; Beaufort et al, 201 1). We also investigated 
the effect of CO2 on nitrate and phosphate utilization by analyz- 
ing the activity of two enzymes involved in nutrient assimilation: 
alkaline phosphatase (APase), and nitrate reductase (NRase). 
Investigating the response of different E. huxleyi strains, and 
potential discrepancies between them, to environmental change 
is central to model the contribution of this ecologically important 
species to the global carbon cycle. 

MATERIALS AND METHODS 
CULTURE CONDITIONS 

Experiments were conducted in diluted batch cultures of E. hux- 
leyi (Lohmann) W. W. Hay and H. P. Mohler, strain NZEH 
(CAWPO 6), isolated in 1992 in the South Pacific Ocean and 
obtained from the Plymouth Culture Collection (UK). Artificial 
sea-water (ASW) was prepared according to Kester et al. (1967) 
with different nitrogen and phosphorus concentrations to achieve 
nutrient (nitrate, phosphatej-replete (R), nitrate-limited (-N) 
and phosphate-limited (-P) conditions (Table 1). Trials were con- 
ducted prior to the experiments to ensure that -N and -P cultures 
reached nutrient limitation at the desired cell density, a den- 
sity below levels that would alter significantly the media carbon 
chemistry. The three nutrient regimes were combined with dif- 
ferent CO2 partial pressures (ranging from 250 to 1314u,atm) 
corresponding to pre-industrial levels and projected values for the 
middle and the end of the century respectively (Table 1). Medium 
carbonate chemistry was adjusted by additions of sodium car- 
bonate (Na2COa) and hydrochloric acid (HC1) to change the 
relative proportion of dissolved inorganic carbon (DIC) species 
and restore total alkalinity (TA) respectively (Riebesell et al, 
2010). The conditions mimicked changes in carbonate chem- 
istry associated with ocean acidification (CO2 increases while 
TA remains constant at ~2268 ± 64.86 u,mol kg -1 ) (Table 1). 
The culture medium was filtered through sterile 0.22 \im 
polycarbonate filters (Millipore® Stericup™ Filter Units). All 
other environmental parameters remained constant through- 
out the experiments: salinity = 34.00 ± 0.40, temperature = 
19.00 ± 0.50°C, 12:12 h lightdark cycle, irradiance = 120.00 ± 
15.00 (imol photons m 2 s _1 under Sylvania Standard F36W/135- 
T8 white fluorescent lighting (Havells Sylvania, Newhaven, UK). 

INCUBATION EXPERIMENTS 

Experiments were conducted in triplicate in 4 L Nalgene® poly- 
carbonate bottles. After the cell inoculation at an initial density 



of 100 cells ml -1 , the bottles were completely filled to minimize 
headspace, closed and sealed with Parafilm® until harvested. A 
blank control bottle (containing no cells) was incubated alongside 
each treatment. A fourth replicate bottle (seeded with the same 
original stock culture and at the same concentration as that used 
in the triplicate bottle experiments) was used for daily monitoring 
of cell density, temperature, pH and irradiance, to avoid opening 
any of the triplicate bottles during the course of the experiment. 
Cell densities at the time of harvest in R cultures were 75,988 ± 
13,159 cells ml -1 depending on the treatments. The R cultures 
were harvested during exponential growth phase and did not 
experience nutrient limitation over the course of the experiment. 
The -N and -P cultures were harvested two days after exponential 
growth stopped (assessed by daily cell counts from the test bot- 
tle). This allowed cultures to be in growth-limiting conditions for 
2 days (cell densities at the time of harvest were 71,587 ± 9250 and 
43,288 ± 14,651 cells ml -1 for -N and -P cultures respectively). 
All cultures were allowed to grow for 8-10 generations, corre- 
sponding to a maximal DIC consumption of 12%. This number 
of generations ensured that almost 100% of the cells in the cul- 
tures experience the study conditions. At time 0 (pre-inoculation) 
and during harvesting (always conducted 3 h after the beginning 
of the light phase), samples were collected from all experimen- 
tal bottles for analysis of carbon chemistry and macronutrient 
concentration in the medium, particulate organic carbon (POC), 
Ca 2+ measurements [for determination of particulate inorganic 
carbon (PIC)], particulate organic nitrogen (PON) and phos- 
phorous (POP), cell density, and scanning electron microscope 
(SEM) imaging. Samples were also collected for APase and NRase 
assays. 

DETERMINATION OF NITRATE REDUCTASE ACTIVITY 

Aliquots of 400 ml were centrifuged (2000 g, 4°C, 15min), and 
the resulting pellets were snap frozen in liquid nitrogen and stored 
at — 80°C. NRase was extracted by adding 500 u,l of a solution 
containing 0.20 M phosphate buffer (pH = 8.20), 1 mM dithio- 
threitol (DTT) and 0.50 M methylenediaminetetra-acetic acid 
(EDTA) to each pellet. The resuspended material was sonicated 
on ice for nine 10-s bursts (30 s intervals between bursts) using 
a VC300 Vibracell sonicator (Sonics and Materials, USA) with a 
20-kHz frequency, 50% duty cycle and an output of 3 (90 W). 
The final extract was centrifuged again (750 g, 4°C, 5 min), and 
the supernatant was used for the enzyme activity determination. 
NRase assays were developed according to Rigobello-Masini et al. 
(2006). Tests were carried out in triplicate in 1 ml at 19°C. The 
reaction mixture contained 100 |xl of crude extract, 10 mM KNO3 
and 2 mM MgS04 and was initiated by the addition of reduced 
nicotinamide adenine dinucleotide (NADH) substrate to a final 
concentration of 0.40 mM. The NRase reaction was stopped after 
15 min with 250 jxl of absolute ethanol at 0°C and with 50 mM 
ZnS04. Activity was estimated based on the final nitrate con- 
centration, indicated by the formation of a red AZO product 
after the simultaneous addition of 100 \il (0.10% weight in vol- 
ume, w/v) sulphanilamide and 100 [xl (0.10% w/v) N-l-naphtyl 
ethylenediamine dihydrochloride (Nicholas and Nason, 1957). 
After these additions, the reaction mixture was centrifuged again 
(21,000 g, 5 min), and the supernatant taken for colorimetric 
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Table 1 | Nutrient concentrations and carbon chemistry parameters of the media at the beginning and at the end of the experiment. 

Nutrient [NO3-] [PO43-] DIC ((imol TA ((jimol pH tota | pC0 2 [HCO3-] [CO3-] C0 2 S2-Ca 

condition (n.molkg _1 ) ((imol kg -1 ) kg - ) kg - ) atm) (jimol kg - ) (|imol kg - ) (iimol kg - ) 



Initial 3 


R 


161.2 


3.5 


1906 


2230 


8.18 


258 


1674 


222 


8 


5.3 


End b 




149.0 (10.4) 


3.0 (0.1) 


1808 (5) 


2141 (4) 


8.21 (0.00) 


225 (2) 


1575 (5) 


225 (1) 


8(0) 


5.4 (0.0) 


Initial 


R 


154.4 


4.1 


2040 


2221 


7.90 


555 


1890 


132 


19 


3.2 


End 




147.6 (0.9) 


2.6 (0.0) 


1937 (65) 


2113 (7) 


7.91 (0.01) 


519(9) 


1793 (5) 


127 (2) 


17 (0) 


3.0 (0.0) 




p> 


156 2 


3 4 


2238 


2330 


7 67 


1073 


2117 


86 


36 


2 1 


End 




14Q 4 to q\ 


1 7inii 


910? (141 


2179 (12) 


7 R4 rn m *i 


1080 (?7\ 


1 ..J .3'. J \ 1 


75 (1 ) 


3fi (1 \ 


1 8 (0 0) 


Initial 


-N 


3.7 


3.6 


1892 


2221 


8.19 


250 


1658 


225 


8 


5.4 


End 




0.3 (0.1) 


2.9 (0.4) 


1734 (4) 


2021 (3) 


8.16 (0.01) 


251 (5) 


1534 (5) 


192 (3) 


8(0) 


4.6 (0.1) 


Initial 


-N 


3.1 


3.5 


2085 


2308 


7.99 


464 


1909 


161 


16 


3.9 


End 




0.1 (0.0) 


2.9 (0.0) 


1917 (4) 


2111 (5) 


7.95 (0.00) 


465 (4) 


1765 (3) 


137 (1) 


16(0) 


3.3 (0.0) 


Initial 


-N 


3.3 


3.5 


2195 


2264 


7.60 


1229 


2082 


72 


41 


1.7 


End 




0.1 (0.0) 


1.6 (0.00) 


1958 (30) 


1994(1) 


7.51 (0.01) 


1358 (25) 


1860 (3) 


52 (1) 


45(1) 


1.3 (0.0) 


Initial 


-P 


158.8 


0.2 


1897 


2217 


8.18 


256 


1667 


222 


9 


5.3 


End 




156.7 (0.21) 


0.0 (0.0) 


1790 (7) 


2097 (8) 


8.18 (0.00) 


243 (3) 


1573 (7) 


209 (2) 


8(0) 


5.0 (0.0) 


Initial 


-P 


159.5 


0.2 


2169 


2363 


7.92 


561 


2004 


147 


19 


3.5 


End 




153.3 (1.5) 


0.0 (0.0) 


1890 (3) 


2109 (8) 


8.01 (0.01) 


394(14) 


1723 (5) 


154 (5) 


13 (1) 


3.7 (0.1) 


Initial 


-P 


152.7 


0.2 


2186 


2241 


7.57 


1314 


2075 


67 


44 


1.6 


End 




151.5 (0.2) 


0.0 (0.0) 


2098 (3) 


2153 (3.3) 


0.57 (0.00) 


1253 (8) 


1992 (3) 


64.8 (0.4) 


42 (0) 


1.6 (0.0) 



"Average blank values at the beginning of the experiment, before the inoculation of the cells. 
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R, nutrient replete; -N, nitrate limited; -P, phosphate limited conditions. 



analysis at 543 nm. Absorbance values were converted to nitrate 
concentration using a calibration curve. Enzymatic activities were 
expressed in enzymatic units per total protein, where one unit of 
the enzyme activity (UEA) catalyzes the conversion of 1 (jimol of 
nitrate to nitrite. Total protein determinations were performed 
with a commercial kit (Thermo Scientific Pierce® BCA Protein 
Assay Kit). 

DETERMINATION OF ALKALINE PHOSPHATASE ACTIVITY 

APase is expressed on the cell surface of E. huxleyi (Dyhrman 
et al., 2006), allowing the activity assay to be performed on whole 
cells. Aliquots (40 ml) of sample were centrifuged at 2000 g and 
19°C for 15min, and the pellets were resuspended in 2.90 ml 
of a solution containing 0.01 M Tris buffer (pH = 9.00), 0.05 M 
MgCl 2 and 0.01 M CaCl 2 . After resuspension, 100 |xl of 13.50 M 
p-nitrophenilphosphate (p-NPP; Sigma) substrate were added, 
and the mixture was incubated at 19°C for 20 min. The reac- 
tion was stopped by the addition of 0.60 ml of 1 M NaOH, and 
samples were centrifuged again (3000 g, 19°C, 15 min) before 
the absorbance of the supernatant at 410 nm was measured. 
Absorbance values were transformed to p-nitrophenol (p-NP) 
concentration using a suitable calibration curve. One unit of 
enzymatic activity corresponds to 1 nmol of p-NP produced per 
10 6 cells min" 1 . 

GROWTH RATE AND C0CC0SPHERE VOLUME 

Growth rate was determined with a standard exponential growth 
equation (Reynolds, 1984): 

tx= (Ln(N t )-Ln(N 0 ))/t (1) 



where Nq and N t are the cell densities at the start and at the 
harvest day respectively, and t corresponds to the length of 
incubation (in days). Cell density and estimated coccosphere 
(cell + coccoliths) volume were determined in triplicates using 
a Beckman Coulter Multisizer III with a 70 u,m aperture. 

C0CC0LITH LENGTH 

Coccolith length was measured from SEM images. For SEM sam- 
ple collection, a 25 mm MF 300 filter was soaked with a drop of 
dilute ammonium hydroxide, and a 0.22 [im polycarbonate fil- 
ter was placed on top. A few drops of culture were placed on 
the polycarbonate filter, and samples were dried on an open Petri 
dish (37°C for 24 h). A section of the top filter was cut out and 
sputter-coated in a Hummer VI-A gold coater, and a grid of 
100 images at 5000 x magnification was taken at a random loca- 
tion on each filter using a LEO 1450 VP SEM with SmartSEM 
V05-1 software. At least 60 coccoliths were measured on con- 
secutive images along their longest axis (defined as coccolith 
length). 

PARTICULATE MATTER ANALYSES 

POC and PON concentrations were measured using a Thermo 
Finnigan Flash EA1112 elemental analyzer with acetanilide stan- 
dards at Plymouth Marine Laboratory (PML). Aliquots of 200 ml 
were filtered through two pre-combusted (400°C, 4h) MF 300 
filters (25 mm glass microfiber 0.70 [im pore size, Fisherbrand). 
Filters were kept at — 20°C until required for analysis and 
fumed with sulphurous acid for 24 h in a desiccator chamber 
to remove inorganic carbon (Verardo et al., 1990). The filters 
were then dried at 60° C for 16 h and pelleted in pre-combusted 
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aluminium foil (EMA; 100 x 30 mm circles) following Hilton 
etal. (1986). 

For PIC analysis, 200 ml of medium was filtered through 
0.20 u,m 47 mm diameter Nuclepore polycarbonate filters, pre- 
viously rinsed twice with 5 ml of dilute ammonium hydroxide 
solution (pH ~9), and washed again three times after filtering. 
Filters were stored in 50 ml Falcon tubes at — 20° C until analy- 
sis. Samples (including blanks) were then weighed, and 15 ml of 
0.10 M nitric acid were added and re-weighed to determine the 
acid volume. The filters were left in acid for 2-3 h with contin- 
uous shaking, after which 500 u,l of the acid leach was removed 
and centrifuged (6500 g for 6 min). A 250 |xl aliquot of the super- 
natant was taken to determine elemental concentrations in a 
Varian Vista Pro ICP-OES. The Ca 2+ per coccolithophore was 
calculated and extrapolated to PIC, assuming that all Ca 2+ on 
the filters originated in CaCC>3 (Fagerbakke et al, 1994). The 
precision of the method, assessed from periodical measurements 
(n = 11) of standards was 1.17% RSD. 

POP was measured using a wet-oxidation method, as 
described by Raimbault et al. (1999). Medium aliquots (200 ml) 
were filtered through a single pre-combusted (400° C, 4 h) MF 300 
filter (25 mm glass microfiber 0.70 [im pore size, Fisherbrand). 
Samples were digested with sodium tetraborate and potassium 
persulphate and autoclaved before analysis in a Segmented Flow 
Auto Analyser (SEAL QuAAtro) at the National Oceanography 
Centre Southampton (NOCS), UK. 

MEDIA CHEMISTRY 

Samples (20 ml) of media were collected for nutrient measure- 
ments by filtration through a 0.22 [im Millex filter (Millipore, 
Billerica, MA, USA) and stored at -20° C until analysis. 
Macronutrient concentrations were determined colorimetrically 
following Hansen and Koroleff (1999) using a Segmented Flow 
Auto Analyser (QuAAtro, SEAL Analytical) at the NOCS (UK). 

Samples for carbonate chemistry were collected in 300 ml 
borosilicate bottles and preserved in the dark with HC1 at a final 
concentration of 2.5 10~ 3 M to prevent microbial growth dur- 
ing storage. These samples were later analysed to determine TA 
and DIC using a Verstatile INstrument for the Determination 
of Total inorganic carbon and titration Alkalinity (VINDTA3C) 
at the NOCS. DIC was analysed using a colorimetric titration 
(coulometer 5011, UIC, USA), and TA was determined using 
a semi-closed cell titration (Dickson et al., 2007). The preci- 
sion of the method, assessed daily from repeated measurements 
{n > 5) on the same batch of seawater, was 2.7 ± 1.6|xmol 
Kg" 1 for DIC and 0.78 ± 0.78(jtmol Kg" 1 for TA. The accu- 
racy was controlled against Certified Reference Materials (from 
A. G. Dickslon, Scripps Institution of Oceanography, USA) mea- 
sured at the beginning and end of each day of analysis apply- 
ing a correction factor obtained from the difference between 
the certified and the measured values. The carbonate system 
was calculated from temperature, salinity, DIC, TA and nutri- 
ents using the "C02SYS" macro (Lewis and Wallace, 1998). 
The equilibrium constants were from Mehrbach et al. (1973) 
and refitted by Dickson and Millero (1987). The KS0 4 con- 
stants were from Dickson (1990), and a seawater pH scale 
was used. 



STATISTICAL ANALYSIS 

One way factor AN OVA was conducted using SPSS 17 (SPSS 
Inc., Chicago IL, USA). Linear correlation factors (r 2 value) were 
calculated using Sigma Plot 11.0 version (Systat Software Inc.). 

LIMITATION OF THE EXPERIMENTAL APPROACH 

Our batch culture experimental design precludes a direct quan- 
titative comparison of PON, POP, PIC, and POC quotas (pmol 
per cell" 1 ) between R and -P or -N cultures. Unlike in R cul- 
tures, growth rates in the -N and -P treatments were not constant 
over the course of the experiment (see Langer et al., 2012; Langer 
et al., who used the same experimental approach). In nutrient- 
limited batch cultures, cells experienced an initial exponential 
nutrient- replete phase followed by a nutrient-limited phase when 
cell division rate decreased. Therefore, cellular quotas include 
cellular PIC and POC produced in both exponential and nutrient- 
limited periods. For this reason, these cannot be compared with 
cellular quotas under nutrient-replete conditions, where growth 
is constant and exponential during the experiment. Therefore, in 
this study, R, -P, and -N cultures are treated as separate experi- 
ments and comparisons of cellular PIC and POC quotas can be 
drawn between CO2 treatments but only within the same nutri- 
ent condition. Within each nutrient-limiting condition, the initial 
nutrient concentrations for all the CO2 conditions were identical 
(see Table 1). Given that in -N and -P experiments cells were har- 
vested 2 days after the end of the exponential growth phase, any 
differences in PIC and POC quotas between the CO2 treatments 
(within each nutrient condition) are only the result of different 
CO2 levels. 

RESULTS 

ENZYMATIC ACTIVITY 

NRase activity was detected only under R conditions, and a lin- 
ear increase in its activity [r 2 = 0.79; f(i,7) = 26.4; p = 0.01] was 
observed with increasing CO2 levels (Figure 1). This increase in 
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FIGURE 1 I Enzymatic response of nitrate reductase (NRase) and 
alkaline phosphatase (APase) to different C0 2 and nutrient scenarios 
in Emiliania huxleyi NZEH. 
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NRase activity was accompanied by a simultaneous increase in the 
cellular POC quota [Figure 2; r 2 = 0.74;_F (2 ,7) = 20.5;p = 0.02]. 
APase activity was only detected in the -P cultures and showed 
a maximum rate of 6.25nmoles of p-NP 10 6 cells min -1 at 
561 (Jtatm CO2. APase activity at 256 and 1314 |xatm CO2 was 77 
and 61% lower than at 561 |xatm C0 2 [Figure 1; F( 2 ,6) = 149.0, 
p < 0.001]. Cellular PIC quota showed a strong correlation with 
APase activity [Figure 3; r = 0.94; F (h7) = 59.1; p < 0.001] but 
not with cellular POC (Figure 3). 

CELL ORGANIC AND INORGANIC MATTER QUOTAS 

In R cultures, both cellular PON and POP quotas increased 
significantly with increasing CO2 [Figure 4A; -F(2,6) = 19.28, 
p = 0.002 (PON); F (2 ,6) = 24.69, p = 0.002 (POP)]. CeUular 
PIC and POC quotas were ~70 and 90% higher respectively 
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at the highest CO2 levels (555 and 1073 (iatm) compared 
with cells grown under ~258 (xatm of CO2 [Figure 4B; 



F(2, 6) = 80.37, 



p < 0.001 (PIC); F| 



(2,6) 



10.65, 



0.01 



(POC)]. In addition, cellular POC quotas increased in parallel 
to PON quotas under different CO2 conditions [Figure 2; 
r 2 = 0.99; F ( i, 7) = 3596.4; p < 0.001]. The PIQPOC ratio 
did not change significantly with rising CO2 [Figure 4C; 
_F( 2 ,6) = 1-41, p = 0.315] although a decreasing trend could be 
observed. 

In -N cultures, nitrate in the culture media was almost com- 
pletely depleted after 8-9 generations (Table 1). The cellular 
PON quota was less than half the concentration of the R cul- 
tures [Table 2; _F (4 , 18) = 9.70, p < 0.001]. Cellular PON quo- 
tas did not vary [Figure5A; F( 2 ,6) = 1-22, p = 0.361] but POP 
quotas increased significantly [Figure 5 A; -F(2,6) = 16.33, p = 
0.004] with increasing CO2. Similarly, PIC quotas did no vary 
[Figure 5B; F(2.6) = 0.35, p = 0.72], but POC quotas increased 
significantly with increasing CO2 levels [Figure 5B; Fq, 6) = 6.14, 
p = 0.035]. A significant decreasing pattern in PIQPOC was 
observed in -N cultures [Figure 5C; F(2,6) = 9.07, p = 0.015] 
with increasing CO2. 

In -P cultures, phosphate in the culture media was almost 
completely depleted after 8-9 generations (Table 1). The cellular 
POP quota was less than half the concentration of the R cul- 



tures [Table 2; F, 



(4,18) 



7.39, p = 0.001] 



increased significantly [Figure 6A; F, 



(2,6) 



Cellular PON quotas 
= 9.34, p = 0.014] but 



POP quotas did not vary [Figure 6A; F (2 , 6 ) = 1.59, p = 0.279] 
with increasing CO2 conditions. The trend in cellular PIC quo- 
tas was not uniform and the most pronounced increase (> 100%) 
was observed at 561 (jtatm of CO2 [Figure 6B; F(2,6) = 37.03, 
p < 0.001]. Cellular POC quotas were ~50% higher at the two 
highest CO2 levels (~561 and 1314u,atm) compared with cells 
grown under 256 (xatm CO2 [F(2,6) = 10.742, p = 0.010]. In gen- 
eral, within each nutrient condition, the trend in PIC and POC 
quotas with CO2 varied when expressed as a function of coc- 
cosphere volume (Table 2). PIQPOC showed a decreasing pat- 
tern with increasing CO2 conditions [Figure 6C; F(2,6) = 8.97, 
p = 0.016]. 

C0CC0LITH LENGTH AND C0CC0SPHERE VOLUME 

Under the highest CO2 level used, coccolith length was ~10% 
higher compared to the lowest CO2 condition for all the 
nutrient treatments [Figure 4D; F(2,567) = 31.7, p < 0.001 (R); 
Figure 5D;F (2 ,567) = 14.5, p < 0.001 (-N); Figure 6D; F (2 , 5 67) = 
43.8, p < 0.001 (-P)]. Coccosphere volume increased by ~39% 
over the CO2 range considered under R conditions [Figure 4D; 
F (2 6 ) = 12.8, p = 0.007]. No differences in coccosphere volume 
were observed under -N conditions for all the CO2 levels tested 
[Figure 5D; F(2,6) = 2.1, p = 0.207]. The highest coccosphere 
volume under -P conditions was observed at 561 (xatm of CO2 
[Figure 6D; F(2, 6 ) = 37.5, p < 0.001]. 

DISCUSSION 

NUTRIENT UTILIZATION: ALKALINE PHOSPHATASE AND NITRATE 
REDUCTASE 

This is the first study to investigate E. huxleyi APase and NRase 
activities under varying CO2 levels. APase was only detected in -P 
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cultures since APase activity is typically enhanced by phospho- 
rus limitation. This enzyme allows phytoplankton to overcome 
phosphorus starvation by hydrolysing phosphate from esters in 
the dissolved organic phosphorus pool. APase activity showed a 
clear response to CO2 partial pressure, and its maximum activ- 
ity was found at 561 u,atm CO2 (Figure 1). Very little is known 
about the effect of CO2 in APase activities and phytoplankton 
(Tanaka et al, 2008; Endres et al, 2013), and the only study 
with E. huxleyi focused on the effect of pH (Xu et al., 2006). 
Since both CO2 and pH covaried in this study, the independent 
effect of these two parameters in APase activity cannot be dis- 
cerned. As APase activity is known to be affected by changes 
in pH (Kuenzler and Perras, 1965), it is likely that pH is one, 
although maybe not the only, parameter affecting APase activity 
in this study. 

A non-uniform response of APase activity to CO2 levels 
was observed, increasing from 256u,atm CO2 {pH = 8.18) to 
561 u,atm CO2 (pH = 7.92) but decreasing at the highest CO2 
conditions (1314 u,atm, pH = 7.57) (Figure 1). This decline in 
activity may be the result of a CO2 threshold exceeding the 
pH upper limit for optimum APase activity. In a study on 
the cyanobacterium species Nodularia spumigena, APase activ- 
ity increased with CO2 (Endres et al, 2013). However, direct 
comparison across studies is not always possible because of the 
different CO2 levels applied. For example, the maximum CO2 
conditions used in their study was ~700 u,atm, probably still 
within the optimal pH range for this species. Additionally the 
thresholds of PH/CO2 tolerance cannot be generalized as APase 
optimum pH is species-specific (Kuenzler and Perras, 1965) 
and also dependent on the species biogeography (Yamada and 
Suzumura, 2010). The decline in APase activity at high CO2 
(1314u,atm) suggests that the competitive ability of E. huxleyi 
to acquire phosphorus may be compromised in future more 
"acidic" oceans. However, the ecological implications of phys- 
iological results should be carefully considered given the high 
degree of genetic diversity among E. huxleyi strains (Iglesias- 
Rodriguez et al., 2006), the marked differences between strain 
maximum activities (Xu et al, 2010; Reid et al., 2011), as well as 
synergies with other levels of ecological organization. 

NRase was only found in nutrient replete (R) cultures, 
when nitrate was present in the medium (Figure 1). Previous 
studies indicate a down-regulation of proteins involved in the 
acquisition and assimilation of inorganic nitrogen after nitrate 
depletion in E. huxleyi (Bruhn et al, 2010). Despite the pres- 
ence of nitrate in the medium, NRase was absent when cul- 
tures grew in -P conditions. Since some regulatory mechanisms 
involved in nitrate assimilation include phosphorylation, phos- 
phate limitation may be impairing this process (Beardall et al., 
1998). NRase increased with rising CO2 levels under R condi- 
tions. The significant positive correlation between cellular POC 
and PON with NRase (Figure 2) could be explained by nitro- 
gen and carbon metabolism being tightly coupled, such that 
a decrease in photosynthetic carbon fixation limits nitrogen 
assimilation (Hipkin et al., 1983; Turpin, 1991). It is possi- 
ble that responses in NRase activity to ocean acidification may 
be species- and possibly strain-specific such as those found 
in higher plants (Fonseca et al., 1997; Matt et al, 2001) and 



phytoplankton species (Xia and Gao, 2005; Rigobello-Masini 
et al., 2006). 

PARTICULATE INORGANIC AND ORGANIC CARBON QUOTAS UNDER 
NUTRIENT-REPLETE CONDITIONS 

Both cellular POC and PIC quotas increased at high CO2 under 
nutrient replete (R) conditions (Figure 4B). This increase was 
accompanied by slight increases in coccosphere volume and 
coccolith size (Figure 4D). The observed increase in cellular POC 
quotas at high CO2 concentrations has been observed in previ- 
ous studies (e.g., Riebesell et al., 2000; Iglesias-Rodriguez et al., 
2008; Langer et al, 2009) and suggests that this species could be 
carbon limited in the present ocean. Considering the low affinity 
of RUBISCO for CO2 and the poorly efficient carbon concen- 
trating mechanisms in E. huxleyi (Paasche, 2002), any increase 
in available CO2 would increase the speed of carbon fixation, 
and thereby cellular POC quotas (Barcelos e Ramos et al, 2010). 
The observed increase in cellular PIC quotas with rising CO2 lev- 
els is in accordance with other studies using E. huxleyi NZEH 
(Iglesias-Rodriguez et al., 2008; Shi et al., 2009) but in disagree- 
ment with Hoppe et al. (2011). Discrepancies between studies 
with the same strain have also been found with strains oth- 
ers than NZEH (Langer et al, 2009; Hoppe et al., 2011) and 
with parameters others than PIC such as growth rate. Several 
explanations could reconcile the different outcomes observed 
between studies using the same strain; for example, differences 
in the experimental set up (i.e., carbon manipulation meth- 
ods, number of generations under the experimental conditions, 
differences in temperature or light). The method of pH manip- 
ulation does not seem to be driving the reported differences 
between experiments according to results by Hoppe et al. (2011) 
and Shi et al. (2009), which revealed similar PIC responses 
to increasing CO2 using bubbling with different CO2 partial 
pressures or acid/base addition. Regarding the number of gen- 
erations exposed to the manipulation, Barcelos e Ramos et al. 
(2010) showed a rapid response of E. huxleyi's metabolic prop- 
erties (including PIC quotas) in response to ocean acidification. 
Finally, any differences in laboratorial culturing conditions, e.g., 
temperature and irradiance could give different outcomes. For 
example, results from this study and those by Iglesias-Rodriguez 
et al. (2008) and Shi et al. (2009) cannot be compared with 
those by Hoppe et al. (2011) because the temperature used in 
the latter was 4-5° C higher. Additionally, synergistic effects can 
give different outcomes; for example, trends in cellular PIC can 
change when CO2 levels are combined with different tempera- 
tures (Borchard et al., 2011) or light levels (Zondervan et al., 
2002). 

Different outcomes between studies using the same strain 
could also be explained by genetic differences between the cul- 
tured strains themselves. For example, there is evidence of shifts 
in phenotypic and genomic properties of strains over time under 
continuous culturing in the laboratory (Lakeman et al., 2009). 
Also, the potential for a strain to evolve properties that deviate 
from those of its original phenotype it is known to increase with 
time of exposure under the new growth conditions (Lakeman 
et al, 2009). Thus, comparisons between studies using the "same" 
strain must be conducted with caution. 



Frontiers in Microbiology | Aquatic Microbiology 



June 2013 | Volume 4 | Article 155 | 8 



Rouco et al. 



E. huxleyi under nutrient limitation and ocean acidification 



PARTICULATE INORGANIC AND ORGANIC CARBON QUOTAS UNDER 
NUTRIENT-LIMITED CONDITIONS 

The ocean is a dynamic system and the physiological response 
of an independent species might differ depending on the 
combination of environmental parameters or stressors to which 
they are exposed, including nutrient limitation. Similarly to 
R cultures, the increase in POC quotas with increasing CO2 
levels in both -N and -P cultures (Figures 5B, 6B) suggests 
that, under nutrient limitation, carbon may be rate-limiting 
for photosynthesis under the CO2 conditions commonly found 
in the open ocean. Our results are however, in disagreement 
with those using nitrogen (Sciandra et al., 2003) and phos- 
phorous (Borchard et al, 2011) limiting conditions in a cal- 
cifying E. huxleyi strain, but in accordance with Leonardos 
and Geider (2005), using a non-calcifying E. huxleyi strain 
under high irradiance. These differences probably result from 
strain-specific responses and/or variations in laboratorial con- 
ditions. In contrast to cellular POC quotas, PIC quotas in - 
N and -P cultures showed a different trend to that observed 
under R conditions. Interestingly, in -N cultures, PIC quotas 
remained constant under the different CO2 conditions suggesting 
that nitrogen metabolism may be decoupled from calcification. 
However, under phosphate limitation a decrease in PIC quotas 
was observed under the highest CO2 conditions (1314u,atm). 
Interestingly, APase activity was tightly correlated with PIC quo- 
tas (Figure 3). However, elucidating the mechanistic effect of 
CO2 and phosphate limitation on calcification requires fur- 
ther work. 

Different cellular elemental stoichiometry was found in the 
E. huxleyi NZEH strain under different environmental scenar- 
ios (Table 2). Stoichiometric mechanisms are known to play 
a very important role in defining the structure of the food 
web in aquatic ecosystems (Elser et al, 2000). Specifically, 
the biochemical composition of phytoplankton is known to 
affect grazing preference (Jones et al, 2002). Thus, changes 
in stoichiometry could influence the grazing-selection pres- 
sure, ultimately determining the prevalence of some strains 
versus others in future oceans. In addition, changes in cellu- 
lar elemental stoichiometry should be taken into consideration 
when predicting the role of E. huxleyi in future biogeochemical 
cycles. 

Unlike in the R manipulations, the PIC:POC ratios decreased 
with CO2 under both -N and -P conditions in the NZEH strain 
(Figures 5C, 6C), comparable to other E. huxleyi strains under 
nitrogen (Sciandra et al., 2003; Mtiller et al., 2012) and phos- 
phorus (Borchard et al, 2011) limitation. This finding is impor- 
tant considering the geographic extent of oligotrophic oceanic 
waters and the importance of PIQPOC ratio in determining 
aggregate formation processes, properties and sinking velocities 
(Armstrong et al., 2002). The data presented here suggest that a 
N- and P-limited population in a high CO2 ocean would have 
a reduced PIQPOC ratio, which would increase the removal of 
CO2 as a combined effect of calcification (a source of CO2) 
and photosynthesis (a sink of CO2) (Frankignoulle et al., 1994). 
However, this could also reduce the net CO2 export to the deep 



ocean associated with the role of coccolith CaC03 in forming 
aggregate as ballast (Armstrong et al., 2002; Biermann and Engel, 
2010). 

COCCOLITH LENGTH AND C0CC0SPHERE VOLUME 

Although coccolith length was correlated with CO2, cocco- 
sphere volume was better correlated with nutrient availability 
(Figures 4D, 5D, 6D). Coccosphere volume increased with ris- 
ing CO2 under R conditions, but an 11% decrease in coc- 
cosphere volume was observed under the highest CO2 level 
in -P cultures, similar to previous studies with the E. huxleyi 
strain PML-B92/11 (Borchard et al, 2011). It is interesting to 
note that, unlike R and -N cultures, changes in POC quotas 
were not associated with changes in coccosphere volume in -P 
cultures. Additionally, maximum PIC quotas were observed at 
561 iiatm CO2 while coccolith length showed the highest values 
at 1314|xatm CO2. Similarly, in -N cultures, increased coccol- 
ith sizes did not correlate with PIC quotas, which were constant 
at all CO2 conditions. These results suggest that under nutrient 
limitation and at the highest CO2 condition, E. huxleyi might 
hold coccoliths bearing less calcite than under the lowest CO2 
levels. 

IMPLICATIONS 

The E. huxleyi strain NZEH, whose coccolith production seemed 
to be resilient to ocean acidification, presents a different response 
to increasing CO2 depending on the nutrient condition. However, 
and similar to what was found under R conditions, this response 
still seems to be different to that observed in other E. hux- 
leyi strains under -N and -P conditions. For the past few years, 
the ocean acidification community is urging the need for mul- 
tiparametric experiments in order to gain a better insight into 
more realistic species-specific responses to environmental pres- 
sure. However, more strain-specific studies are also necessary in 
order to predict and understand the direction of future changes 
with a degree of certainty. This information is important to 
improve parameterizations in diagnostic and prognostic of global 
biogeochemical models. 
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